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CN I Abstract: We suggest an alternative mathematical model for the electron in 

dimension 1-1-2. We think of our (l-|-2)-dimensional spacetime as an elastic con- 
' fH '' tinuum whose material points can experience no displacements, only rotations. 

Q^, This framework is a special case of the Cosserat theory of elasticity. Rotations 

of material points are described mathematically by attaching to each geomet- 
ric point an orthonormal basis which gives a field of orthonormal bases called 
C^ . the coframe. As the dynamical variables (unknowns) of our theory we choose a 

coframe and a density. We then add an extra (third) spatial dimension, extend 
our coframe and density into dimension l-t-3, choose a conformally invariant La- 
grangian proportional to axial torsion squared, roll up the extra dimension into 
a circle so as to incorporate mass and return to our original (l+2)-dimensional 
spacetime by separating out the extra coordinate. The main result of our pa- 
QQ ' per is the theorem stating that our model is equivalent to the Dirac equation 

^^ , in dimension 1-f 2. In the process of analyzing our model we also establish an 

^T) • abstract result, identifying a class of nonlinear second order partial differential 

f * I equations which reduce to pairs of linear first order equations. 

o 

^^ ' 1. Introduction 

In this paper we consider an electron living in (l+2)-dimensional Minkowski 
■ r-n , spacetime M^+^ with coordinates a;", a = 0, 1, 2, and metric gais — diag(— 1, -1-1, +1). 

^\f ' The reduction of spatial dimension from 3 to 2 makes dealing with spin easier as 

M , we have only two possibilities, spin up and spin down. At a technical level this 

reduction of spatial dimension manifests itself in the fact that we do not need a 
4-component complex bispinor for describing the electron, just a 2-component 
complex spinor. 

The Dirac equation in M^+'^ is 

K^b(za + A)„±ma\b]r;'' = 0. (1.1) 
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Here m is the electron mass, a" are Pauli matrices (see (IA.2I) . (jA.Sp ). da = d/dx'^ 
and Aa is a given external real electromagnetic field. The tensor summation in- 
dex a runs through the values 0,1,2, the spinor summation index b runs through 
the values 1,2 and the free spinor index d runs through the values 1,2. The 
spinor field r] : M^+^ —^ C^ is the dynamical variable (unknown quantity). The 
two choices of sign give two versions of the Dirac equation corresponding to spin 
up and down. 

Equations (jl.ip are, of course, a special case of the Dirac equation in dimen- 
sion 14-3. The latter is a system of four complex equations for four complex 
unknowns and if one looks for solutions which do not depend on x^ then this 
system splits into a pair of systems (|l.ip . 

Throughout this paper all fields are assumed to be infinitely smooth with 
no assumptions on their behavior at infinity. We focus on understanding the 
geometric meaning of equation (jl.ip rather than on fitting it into the framework 
of operator theory. 

We suggest a new geometric interpretation of equation (II. ip . The basic idea is 
to view our (l+2)-dimensional spacetime as an elastic continuum whose material 
points can experience no displacements, only rotations, with rotations of different 
material points being totally independent. The idea of rotating material points 
may seem exotic, however it has long been accepted in continuum mechanics 
within the Cosserat theory of elasticity [7] . This idea also lies at the heart of the 
theory of teleparallelism (= absolute parallelism = fernparallelismus), a subject 
promoted by A. Einstein and E. Cartan [51 [mil2j . With regards to the latter it 
is interesting that Cartan acknowledged [3] that he drew inspiration from the 
'beautiful' work of the Cosserat brothers. 

An elastic continuum with no displacements, only rotations, is, of course, a 
limit case of Cosserat elasticity. The other limit case is classical elasticity with 
displacements only and no (micro)rotations. 

Rotations of material points of the (l-|-2)-dimensional elastic continuum are 
described mathematically by attaching to each geometric point of Minkowski 
spacetime 1^-"^+^ an orthonormal basis, which gives a field of orthonormal bases 
called the frame or coframe, depending on whether one prefers dealing with 
vectors or covectors. Our model will be built on the basis of exterior calculus so 
for us it will be more natural to use the coframe. 

Our model is described in Section [5] Subsequent sections contain mathemati- 
cal analysis culminating in Theorem [T] (see Section [S]) which establishes that our 
model is equivalent to the Dirac equation (|l.ip . 

The mathematical model presented in Section [2] is quite simple. However, 
seeing that this model generates the Dirac equation (|l.ip is not easy. The main 
difficulties are as follows. 

— The dynamical variables in our model and the Dirac model are different. We 
will overcome this difficulty by performing a nonlinear change of dynamical 
variables given by the explicit formulas p.2p - p.4p . 

— We incorporate mass and electromagnetic field into our model by means 
of a Kaluza-Klein extension, i.e. by adding an extra spatial dimension and 
then separating out the extra coordinate x^. Now, our field equation (Euler- 
Lagrange equation) will turn out to be nonlinear so the fact that it admits 
separation of variables is nontrivial. We will establish separation of variables 
by performing explicit calculations. We suspect that the underlying group- 
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theoretic reason for our nonlinear field equation admitting separation of vari- 
ables is the fact that our model is U(l)-invariant, i.e. it is invariant under the 
multiplication of the spinor field by a complex constant of modulus 1. Hence, 
it is feasible that one could perform the separation of variables arguments 
without writing down the explicit form of the field equation. 
— Our field equation will be second order so it is unclear how it can be reduced 
to a first order equation (jl.ip . This issue will be addressed in Appendix [BJ 
Namely, in this appendix we prove an abstract lemma showing that a certain 
class of nonlinear second order partial differential equations reduces to pairs 
of linear first order equations. To our knowledge, this abstract lemma is a 
new result. 

Our paper is a development of the publication [B] where a similar model was 
suggested for a massless fermion (neutrino). 



2. Our model 

The coframe i? is a triple of orthonormal covector fields t?-', j = 0,1,2, in M"'^"'"^. 
Each covector field d^ can be written more explicitly as "d^ a where the tensor 
index a = 0, 1, 2 enumerates the components. Of course, orthonormality is un- 
derstood in the Lorentzian sense: the inner product ^^ • ^'^ = g"'^ 'd^ a'&^ p is —1 
if j = fc = 0, -1-1 if j = A; = 1 or J = fc = 2, and zero otherwise. 

The orthonormality condition for the coframe can be represented as a single 
tensor identity 

g = o,k^^®d'' (2.1) 

where 

ojfc = o^*^ := diag(-l, +1, +1). (2.2) 

For the sake of clarity we repeat formula (j2.ll) giving tensor indices explicitly 
and performing summation over frame indices explicitly: gap = —'d'^a'&^ p + 
'&^a''9^p + 'd'^a'&'^p where a and /3 run through the values 0, 1,2. We view the 
identity (12. 1[) as a kinematic constraint: the covector fields d^ are chosen so that 
they satisfy (j2.ip . which leaves us with three real degrees of freedom at every 
point of M^+^. If one views d^ a as a 3 x 3 real matrix- function, then condition 
(|2.ip means that this matrix-function is pseudo-orthogonal, i.e. orthogonal with 
respect to the Lorentzian inner product. 

We choose to work with coframes satisfying conditions 

Aet§\ = +l>Q, i9%>0 (2.3) 

which single out coframes that can be obtained from the trivial (aligned with 
coordinate lines) coframe "d^ a — 5^ a by proper Lorentz transformations. 

As dynamical variables in our model we choose the coframe 'd and a positive 
density p. Our coframe and density are functions of coordinates a;", a = 0,1,2, 
in M^+^. At a physical level, making the density p a dynamical variable means 
that we view our continuum more like a fluid rather than a solid: we allow 
the material to redistribute itself so that it finds its equilibrium distribution. 
Note that the total number of real dynamical degrees of freedom contained in 
the coframe i) and positive density p is four, exactly as in a two-component 
complex- valued spinor field rj. 
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In order to incorporate into our model mass and electromagnetic field we 
perform a Kaluza-Klein extension: we extend our original (l+2)-dimensional 
Minkowski spacetime M^+^ to (l+3)-dimensional Minkowski spacetime M^+^ 
by adding the extra spatial coordinate x^ . The metric on M^+'^ is ga/3 — 
diag(— 1, +1, +1, +1). Here and further on we use bold type for extended quan- 
tities. Say, the use of bold type in the tensor indices of ga/3 indicates that a. and 
/3 run through the values 0,1,2,3. 

We extend our coframe as 

^'o. = {^'A , J = 0, 1, 2, ^\ = (^^\ (2.4) 

where the bold tensor index a. runs through the values 0,1,2,3, whereas its 
non-bold counterpart a runs through the values 0, 1, 2. In particular, the Oq. in 
formula (J2.4I) stands for a column of three zeros. 

Our original (l+2)-dimensional coframe 1?, which was initially a function of 
(a;", x^, x^) only, is now allowed to depend on x^ in an arbitrary way, as long as 
the kinematic constraint (j2.ip is maintained. Our only restriction on the choice 
of extended (l-t-3)-dimensional coframe 'd is the condition that the last element 
of the coframe is prescribed as the conormal to the original Minkowski spacetime 
Mi+2, see last formula ((2^ . 

We also extend our positive density p allowing arbitrary dependence on x^ . 
We retain the non-bold type for the extended p. 

The coframe elements i9J are different at different points x G M^^'' and this 
causes deformations. As a measure of these "rotational deformations" we choose 
axial torsion which is the 3-form defined by the formula 

T^'^ := iojk-iJJ A d-d^ (2.5) 

where Ojk = oJ'^ := diag(— 1, +1, +1, +1) (compare with formula (j2.2p ) and d 
denotes the exterior derivative on M^+'^. Here "torsion" stands for "torsion of the 
teleparallel connection" with "teleparallel connection" defined by the condition 
that the covariant derivative of each coframe element iP is zero; see Appendix A 
of [2 for a concise exposition. "Axial torsion" is the totally antisymmetric part 
of the torsion tensor (j8.2p . 

We choose the basic Lagrangian density of our mathematical model as 

i(,9,p):=||T-||V (2.6) 

where HT^'^jp = ■jfT^'^^^T^^^g^'^g'^-^g'^'^. The main motivation behind the 
choice of Lagrangian density (j2.6p is the fact that it is conformally invariant: 
it does not change if we rescale the coframe as i9J H> e^^i^J, metric as ga/g 1— >■ 
e^'^ga/3 and density as, p ^^ e^^p where h : M^+^ ^ M is an arbitrary scalar 
function. At this point it is important to note that out that our Kaluza-Klein 
extension procedure does not actually allow for conformal rescalings because the 
last formula (12. 4p is very specific. Thus, our logic is that we choose a Lagrangian 
density (12. 6p which would be conformally invariant if not for the prescriptive 
nature of the Kaluza-Klein construction. This is in line with the view that mass 
breaks conformal invariance. The electron mass m will appear below in formulas 

(I^jtd and isrm . 
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Substituting (1^ into (^3]) we get 

rpax ^ j.ax _ ^3 ^ jj^^ ^2.7) 

where 

T^'^ :== \ojk^^ A dz?*^ (2.8) 

o 

is the axial torsion in original (l+2)-diniensional spacetime (with d now denot- 
ing the exterior derivative on M^+^) and D^d is the 2-form 

D^-d := io.fcz?^" A 93^^ (2.9) 

The 2-form D^-d characterizes the rotation of the coframe "d as we move along 
the coordinate x^ and is, in effect, an analogue of angular velocity. 

Substituting (12. 7p into (|2.6p we rewrite our basic Lagrangian density as 

L{d,p):={\\T'^^^ + \\D,df)p. (2.10) 

We now incorporate the electron mass m into our model by imposing the 
periodicity conditions 

d{x°,x^,x^,x^ +TT/m) =i9(x",x\x2,x^), (2.11) 

p(a;^a;^a;^x^^-7^/TO) = p(x",x\x^ a;^). (2.12) 

Conditions (12. lip and (j2.12p mean that we make the coordinate x^ cyclic with 
period — . In other words, we effectively roll up our third spatial dimension into 
a circle of radius -^ . 

Ira 

Finally, we incorporate the prescribed electromagnetic (co)vector potential A 
into our model by formally mixing up the partial derivatives appearing in the 
definition of axial torsion (12.81) as 



9„^aa+m-Xa3, a = 0,1, 2. (2.13) 

As a result, our Lagrangian density (j2.10p turns into 

L(^,p):=(||rrf + p3^f)p, (2.14) 

where 

T^^ --T^^^rn-^AKD^d. (2.15) 

Let us summarize the above construction. The Lagrangian density that we 
shall be studying is given by formula (J2.14I) where the 3-form T^ and 2-form 
D-^-Q are defined by formulas (|2.8p . (|2.9p and (|2.15p . The corresponding action 
(variational functional) is 

S{-d,p):^f L{i3,p)dx°dx^dx^dx^; (2.16) 

JMI + 3 

of course, the integral in (j2.16p need not converge as we will be using it only 
for the purpose of deriving field equations (Euler-Lagrange equations). Our dy- 
namical variables are the coframe i? and density p which live in the original 
(l-|-2)-dimensional spacetime but depend on the extra spatial coordinate x^. We 
seek solutions which are periodic in x^, see formulas (|2.11l) and (|2.12p . 
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Our field equations are obtained by varying the action (I2.16P with respect to 
the coframe i!) and density p. Varying with respect to the density p is easy: this 
gives the field equation ||T|'^|p + ||£)3i?|p = which is equivalent to L{-d, p) = 0. 
Varying with respect to the coframe d is more difficult because we have to 
maintain the kinematic constraint (|2.ip . A technique for varying the coframe 
with kinematic constraint (I2.ip was described in Appendix B of [3] but we do 
not use it in the current paper. 



3. Switching to the language of spinors 

As pointed out in the previous section, varying the coframe subject to the kine- 
matic constraint (|2.ip is not an easy task. This technical difficulty can be over- 
come by switching to a different dynamical variable. Namely, it is known that 
in dimension 1-1-2 a coframe i!) and a positive density p are equivalent to a 2- 

component complex- valued spinor field ^ = ^^ = ^2 I satisfying the inequality 

e'73abe > 0. (3.1) 

The explicit formulas establishing this equivalence are 

P = CVsafcC', (3.2) 

^"a = p-'C<Jaabe, (3.3) 

(^1 + i^2)^ ^ p-h''a^i^^e<^o.cd^''- (3.4) 

Here a are Pauli matrices and e is "metric spinor" (see (jA.ip - (jA.3|) ). the free 
tensor index a runs through the values 0, 1, 2, and the spinor summation indices 
run through the values 1,2 or 1,2. The advantage of switching to a spinor field 
^ is that there are no kinematic constraints on its components, so the derivation 
of field equations becomes straightforward. 

Formulas p.2l) - p.4l) are a variant of those from [B^: in [B^ these formulas 
were written for dimension 3, i.e. for 3-dimensional Euclidean space, whereas in 
the current paper we write them for dimension 1-1-2, i.e. for (l-|-2)-dimensional 
Minkowski spacetime. Both the formulas from 6 and formulas (j3.2D - (j3.4p are a 
special case of those from [8] . 

Remark 1. The right-hand sides of formulas p.2p - p.4p are invariant under the 
change of sign of ^. Hence, the correspondence between coframe and positive 
density on the one hand and spinor field satisfying condition (13. ip on the other 
is one to two. A spinor field is, effectively, a square root of a coframe and a 
density. The fact that the spinor field has indeterminate sign does not cause 
problems as long as we work on a simply connected open se1p, such as the whole 
Minkowski space M^+^. Note that a similar issue (extraction of a single- valued 
"square root" of a tensor) arises in the mathematical theory of liquid crystals [5] . 



^ Here and further on the notions of openness and connectedness of subsets of M^+^ 
understood in the Euclidean sense, i.e. in terms of a positive 3-dimonsional metric. 
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We now need to express the differential forms (J2.8L (j2.9p and (j2.15p via the 
spinor field f . This is done by direct substitution of formulas (J3.2l) - (j3.4p giving 



(*D3§)a = -T- -2 , (3.6) 

34 CT3cd4 

The tensor summation index a in formulas p.Sp and p.7p and the free tensor 
index a in formula p.6p run through the values 0, 1, 2. Formulas p.Sp and (13.61) 
are, of course, a variant of those from [5] : we simply turned 3-dimensional Eu- 
clidean space into (l+2)-dimensional Minkowski space and replaced the extra 
coordinate x" with the extra coordinate x^. 

Substituting formulas (13. 7p and (13. 6p into (I2.14p we arrive at the following 
self-contained explicit spinor representation of our Lagrangian density 

HO ^ 



iiC'^'^abida + m-^A^d3)e - ^''^'^abida + m-^A^ds)^)]'' 

+ Ikf ^adb^sC' - i'laabd^nWl- (3-8) 



Here and further on we write our Lagrangian density and our action as L(S,) 
and S{^) rather than L{d, p) and S{'d, p), thus indicating that we have switched 
to spinors. The spinor field ^ satisfying condition p.ip is the new dynamical 
variable. 

The field equation for our Lagrangian density p.8p is 

- i*D3d)a,a"abd3^'' - a^afc^s ((^^s^)^^')) ^ p'^ La^ab^" - (3.9) 

where the quantities *T^, *D3-d, p and L are expressed via the spinor field £, in 
accordance with formulas (|3.7p . (13.61) . p.2p and p.8p . 

We seek solutions of the field equation (13. 9p which satisfy the periodicity 
condition 

C(x°, x\x'^,x^+ tt/to) = ^{x°,x\x^,x^), (3.10) 

or the antiperiodicity condition 

^ix°,x\x^,x^ + n/m) = -C(x°,x\ a:^ a;^). (3.11) 

The above periodicity/antiperiodicity conditions are our original periodicity con- 
ditions (j2.1ip and (I2.12p rewritten in terms of the spinor field. The splitting into 
periodicity/antiperiodicity occurs because the spinor field corresponding to a 
coframe and a density is determined uniquely modulo sign, see Remark [1] 
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4. Separating out the coordinate x^ 

Our field equation (J3.9I) is higlily nonlinear and one does not expect it to admit 
separation of variables. Nevertheless, we seek solutions of the form 

C(x^ x\ x^ a;3) = 7y(a;0, x\x^) e^™-'. (4.1) 

Note that such solutions automatically satisfy the antiperiodicity condition p. lip : 
the coframe corresponding to a spinor field of the form (j4.ip experiences one full 
turn (clockwise or anticklockwise) in the ('!9^,'!9^)-plane as x^ runs from to — . 
Substituting formula (gj]) into dSJ]), dSS]), ^^ and (jSH) we get 






*^A± = ^-r— — — , (4.2) 



{*D3'd)a = ±—^^ —, 4.3 

P = ta^abv", (4.4) 



^±(.)-- '' 



[[Urf'j"db{id±AUT^''-7fa''ab{idTA)^t)]^ - {vit^zdwff) (4-5) 

where the signs agree with those in (J4.1D (upper sign corresponds to upper sign 
and lower sign corresponds to lower sign). 

Note that the quantities (|4.2p - (|4.5p do not depend on x"^, which simplifies the 
next step: substituting (j4.ip into our field equation (|3.9p and dividing through 
by the common factor e^*™^ we get 

^ ((*Tr±)a"a6(«9 ± A)„ry'' + g'' ab{id ± A)„((*Tr±)ry' 

32m2 
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o\brt ~ p-^L±n:,abv'' = 0. (4.6) 



Observe that formulas (j4.2p - (|4.6p do not contain x^ . Thus, we have shown 
that our field equation (13. 9p admits separation of variables, i.e. one can seek 
solutions of the form (j4.ip . 

Consider now the action 

S±iri):= [ L±{Tj)dx°dx^dx^ (4.7) 

where L±{'i]) is the Lagrangian density (|4.5p . It is easy to see that equation (|4.6p 
is the field equation (Euler-Lagrange equation) for the action (14. 7p . 

In the remainder of the paper we do not use the explicit form of the field 
equation (J4.6L dealing only with the Lagrangian density (j4.5p and action (j4.7p . 
We needed the explicit form of field equations, (13. 9p and (|4.6p . only to justify 
separation of variables. 
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We give for reference a more compact representation of our Lagrangian den- 
sity (14. Sp in terms of axial torsion T^ (see formula (j4.2[) ) and density p (see 
formula ((i^ ): 

i±(^) = -((*rr±)'-ym^)p. (4.8) 

Of course, formula (|4.8p is our original formula (|2.14p with a;'^ separated out. 
The choice of dynamical variables in the Lagrangian density (|4.8I) is up to the 
user: one can either use the x^-independent spinor field rj or, equivalently, the 
corresponding a:'^-independent coframe and a;'^-independent density (the latter 
are related to rj by formulas (I3.2p - (l3.4p with ^ replaced by rf). 



5. Main result 

Let Drs be the linear differential operator mapping undotted spinor fields into 
dotted spinor fields in accordance with formula 

Tj H^ DrsV ^ o-"db{ida+rAa)ri'' + sma^abv'' (5-1) 

where the tensor summation index a runs through the values 0,1,2 and the 
letters r and s take, independently, symbolic values ± (as in Drs) or numerical 
values ±1 (as in the RHS of formula (|5.1|) ). depending on the context. 
The main result of our paper is 

Theorem 1. Let Q be an open subset o/ M^+^ and let rj : f2 ^ C'^ be a spinor 
field satisfying the condition 

t^^iabrf > (5.2) 

(compare with iS.l]) ). Then r/ is a solution of the field equation for the Lagrangian 
density L+ if and only if it is a solution of the Dirac equation D^-^^rj = or the 

Dirac equation D^ rj = 0, and a solution of the field equation for the Lagrangian 

density L_ if and only if it is a solution of the Dirac equation D ^.r] = or the 

Dirac equation D rj = 0. 

Proof. Put 

ir..(^) := \[fi^<j"ab{tda,+rA^)rj'~rj''a''ab{ido.-rA^)7f] +smta\iif. (5.3) 

This is the Lagrangian density for the Dirac equation DrsTj = 0. Formula (|5.3p 
can be rewritten in more compact form as 

Lrs{^)^(-\*TTr + sm)p (5.4) 

where *T|^, r = ±, is the Hodge dual of axial torsion defined by formula (14.21) 
and p is the density defined by formula (|4.4I) . Comparing formulas (|4.8p and 
(|Ei)) we get 

T ( ^ 32m Lr+irj)Lr-{r]) 

9 Lr+[r]) - Lr-[ri) 

Note that the denominator in the above formula is nonzero because condition 
can be equivalently rewritten as Lr+{ri) > Lr-{ri). 
The result now follows from formula (|5.5p and Lemma [1] (see Appendix [BJ. D 
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6. The sign in the inequality (|3.ip 

In Section [31 when switching to the language of spinors, we chose to work with 
spinor fields ^ satisfying the inequality (13.11) . It is natural to ask the question 
what happens if we choose to work with spinor fields ^ satisfying the inequality 

IVsdbl' < 0. (6.1) 

One can check that in this case all our arguments can be repeated with mi- 
nor changes. Namely, in dimension 1+2 a coframe •& and a positive density p 
are equivalent to a 2-component complex-valued spinor field ^ satisfying the 
inequality (|6.ip . with this equivalence described by a slightly modified version 
of formulas p.2p - p.4p . In the end we get an analogue of Theorem [1] for such 
spinors. 

In fact, there is no need to repeat our arguments because there is a bijection 
between spinor fields ^ satisfying the inequality (13.11) and spinor fields ^ satisfying 
the inequality (|6.1|) : 

e ^ r = e-Vsdbf, i^e = e^Vsdbl". (6.2) 

We do not view the transformation (j6.2p as physically significant because 
the primary dynamical variables in our model are coframe and positive density, 
not the spinor field. We view the spinor field merely as a convenient change of 
dynamical variables. If two different spinor fields correspond to the same coframe 
and positive density we interpret them as the same particle. In group-theoretical 
language this means that our model is built on the basis of the pseudo-orthogonal 
group S0(1, 2) rather than the spin group Spin(l, 2). 



7. Plane wave solutions 

In this section we construct a special class of explicit solutions of the field equa- 
tions for our Lagrangian density (J2.14I) . This construction is presented, initially, 
in the language of spinors and under the additional assumption that the elec- 
tromagnetic covector potential A is zero. 
We seek solutions of the form 

^{x°,x\x^,x^) = e-*(p-^+™^')c (7.1) 

where p = {pq,P2tP3) is a real constant covector, r takes the values ±1 and 
^ 7^ is a constant spinor. We shall call solutions of the type (j7.ip plane wave. 
In seeking plane wave solutions what we are doing is separating out all the 
variables, namely, the original variables x = {x^ , x^ , x'^) (coordinates on M^+^) 
and the extra variable x^ (Kaluza-Klein coordinate). 

As usual, our spinor field ^ is assumed to satisfy the inequality p.ip . As 
explained in Sectional this assumption does not lead to the loss of solutions. 

Our field equation p.9p is highly nonlinear so it is not a priori clear that one 
can seek solutions in the form of plane waves. However, plane wave solutions 
are a special case of solutions of the type (|4.1I) and these have already been 
analyzed in preceding sections. Namely, Theorem [T] gives us an algorithm for the 
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calculation of all plane wave solutions (j7.ip by reducing the problem to Dirac 
equations 

DrsV = (7.2) 

for the x^ -independent spinor field 

rj{x\x\x^)^e-'P-'=C- (7.3) 

Here r is the same as in formula (17. ip . i.e. a number taking the values ±1, and s 
is another number, also taking, independently, the values ±1. By Drs we denote 
the differential operators (15. ip . 

Clearly, a Dirac equation (J7.2I) has a nontrivial plane wave solution rj if and 
only if the momentum p satisfies the condition Upp + m^ — 0, so p is timelike. 
Our model is invariant under proper Lorentz transformations of coordinates 
{x'~' , x-^ , x'^) so without loss of generality we can assume that 

Pi^P2= 0. (7.4) 

Combining formulas ((5?T|) . (|A.2p . (|A.3p . (fTS)) and (fTi)) we see that the Dirac 
equation (|7.2I) takes the form 



-Po + sm \ /C^ 
-po - smj \( 



2,-0. (7.5) 



Equation (j7.3p has a nontrivial solution satisfying the inequality p.ip only if 

Po = sm (7.6) 

with the corresponding C given, up to scaling by a nonzero complex factor, by 
the formula 

C^ = (J) • (7-7) 

Combining formulas (|7.ip . (17.41) . (|7.6p and (|7.7p we conclude that our model 
admits, up to a proper Lorentz transformation of the coordinate system in M^+^ 
and complex scaling, four plane wave solutions and that these plane wave solu- 
tions are given by the explicit formula 

Here the numbers r and s can, independently, take values ±1. 

Let us now rewrite the plane wave solutions (17. 8p in terms of our original 
dynamical variables, coframe ^ and density p. Substituting formulae (|A.2p . (jA.Sp 
and (fTS)) into formulae (l3^ - (IO)) we get p = 1, i?°a = S"a and 

\ / 

cos2m{sx^ + rx^) , 'd'^a = - sin2TO(sa;'' + ra;^) | . (7.9) 

sm2m{sx'^ + rx^) I \ cos2m{sx^ + rx^) 

In order to distinguish the two spins we fix x^ and examine how the covectors 
-d^ and i?^ evolve as a function of time x" . We say that spin is up if the rotation 
is counterclockwise and spin is down if the rotation is clockwise. Examination 
of formula (j7.9p shows that we have spin up if s = -1-1 and spin down if s = —1. 
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We will now establish which of the solutions (I7.9p describe the electron and 
which describe the positron. Let us introduce a weak constant positive electric 
field, < Aq < m and Ai = A2 = 0. Then we can repeat the calculation leading 
up to formula (j7.9p . only now we get 

i?^^ — cos 2[{sm — rAQ)x^ + rmx^] 
y sin2[(sTO — rAQ)x^ + rnix^] 

I 

i?2q = - sin2[{ sm - rAo)x° +rm,x'^] | . (7.10) 

y cos 2[{sm — rylo)a;° + rmx^^ 
We define quantum mechanical energy as 

e -.^ \sm - rAo\ (7.11) 

which is half the angular frequency (as a function of time x'^) of the solu- 
tion (|7.10p . We say that we are dealing with an electron if e < to and with 
a positron ii e > m. Examination of formula (j7.1ip shows that we are looking at 
an electron if the signs of r and s are the same and at a positron if the signs of 
r and s are opposite. This means that the electron is described by a wave trav- 
eling in the negative x'^-direction whereas the positron is described by a wave 
traveling in the positive I'^-direction. 

Our classification of plane wave solutions is summarized in Table [TJ 

Table 1. Classification of solutions 1 17.91 1 







s = +l 


s = -1 


r - 


= -1-1 


Electron with spin up 


Positron with spin down 


r - 


= -1 


Positron with spin up 


Electron with spin down 



8. Discussion 

8.1. Distinguishing the electron from the positron. The mathematical model pre- 
sented in this paper allows us to clearly distinguish the electron from the positron. 
This is achieved by using the coframe and positive density as our primary dy- 
namical variables rather than the more traditional spinor field. In other words, 
as explained in the end of Section |6l our model is built on the basis of the 
pseudo-orthogonal group S0(l,2) rather than the spin group Spin(l,2). 



8.2. Problem of vanishing density. The only technical assumption in our analysis 
is that the density p does not vanish. Rephrased in terms of the spinor field, this 
assumption reads as 

e'J3ate + 0, (8.1) 

compare with p.ip and (|6.ip . We do not know how to drop the assumption (|8.ip . 
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8.3. Curved spacetime. One of the advantages of our mathematical model is 
that it does not use covariant differentiation (only exterior differentiation) so the 
generalization to the case of a curveqj (l+2)-dimensional spacetime is absolutely 
straightforward. Covariant derivatives appear only when we switch from coframe 
and density to a spinor field. All our analysis, including Theorem [U carries over 
to the case of curved spacetime. We chose our (l+2)-dimensional spacetime to 
be flat only to make the exposition clearer. 



8.4. Exclusion of gravity. We assumed the (l+2)-dimensional metric g to be 
prescribed (fixed) and the coframe i9 to be chosen so as to satisfy the kinematic 
constraint (|2.ip . As explained in subsection l8.3l the fact that we chose the metric 
g to be Minkowski is irrelevant and all our analysis carries over to the case of 
an arbitrary Lorentzian metric in dimension 1+2. The important thing is that 
the metric g is not treated as a dynamical variable. This means that we chose 
to exclude gravity from our model. 

On the other hand, in teleparallelism it is traditional to view the metric 
as a dynamical variable. In other words, in teleparallelism it is customary to 
view (j2.ip not as a kinematic constraint but as a definition of the metric and, 
consequently, to vary the coframe ^ without any constraints. This is not surpris- 
ing as most, if not all, authors who contributed to teleparallelism came to the 
subject from General Relativity. 

It appears that the idea of working with a coframe subject to the kinematic 
constraint (12.11) is new. 



8.5. Our choice of Lagrangian. We chose a very particular Lagrangian den- 
sity (j2.6p containing only one irreducible piece of torsion (axial) whereas in 
teleparallelism it is traditional to choose a more general Lagrangian containing 
all three pieces (axial, vector and tensor) of the torsion tensor 

T^'^:=Ojki9J«)di?'', (8.2) 

see formula (26) in 9^. Note that when Einstein introduced teleparallelism [5] 
he neglected the axial piece (|2.5p completely. 

In choosing our particular Lagrangian density (12.61) we were guided by the 
principles of conformal invariance, simplicity and analogy with Maxwell's theory. 
The analogy with Maxwell's theory is that we characterize the field strength by 
a differential form, replacing the electromagnetic tensor (2-form) by axial torsion 
(3- form). It appears that the Lagrangian density (j2.6p was never examined. 



8.6. Density as a dynamical variable. We took the positive density of our con- 
tinuum to be a dynamical variable whereas in teleparallelism the tradition is to 
prescribe it as p = ^J\ deig\ . Taking p to be a dynamical variable is, of course, 
equivalent to introducing an extra real positive scalar field into our model. It 
appears that the idea of making the density a dynamical variable is also new. 



■^ Here "curved" refers to the curvature of the Levi-Civita connection generated by the 
metric g, as is customary in General Relativity. 
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8.1. Electron in dimension 1+3. The major outstanding issue is whether we can 
reformulate the Dirac equation in dimension 1+3 using our approach. This would 
mean starting from (l+3)-dimensional spacetime, performing a Kaluza-Klein ex- 
tension to dimension 1+4, choosing the conformally invariant Lagrangian density 
(|2.6p and so on, as described in Section [21 

It seems that the equation we get starting from (l+3)-dimensional spacetime 
and performing the construction described in Section [2] is not the Dirac equa- 
tion in dimension 1+3. Our analysis is heavily dependent on dimension and, 
when starting from (l+3)-dimensional spacetime, we do not appear to get a 
factorization of the Lagrangian density of the type (j5.5p . 

However, the equation we get in dimension 1+3, although nonlinear, seems 
to be very similar to the Dirac equation. The natural way of testing how close 
our equation is to the Dirac equation would be to calculate the energy spectrum 
of the electron in a given static electromagnetic field, starting with the case of 
the Coulomb potential (hydrogen atom). 



8.8. Similarity with the Ashtekar-Jacobson-Smolin construction. The analysis 
presented in our paper exhibits certain similarities with Ti jlO; in that a 3-dimen- 
sional (or, in our case, (l+2)-dimensional) coframe 1} is used as a dynamical 
variable and that a second order partial differential equation is reduced to a first 
order equation. 



A. Notation 

Our notation follows [SlISlfTB] . The only difference with [51IT5] is that in the latter 
the Lorentzian metric has opposite signature. In [6] the signature is the same as 
in the current paper, i.e. the (l+3)-dimensional metric has signature — \- H — \- . 

We use Greek letters for tensor (holonomic) indices and Latin letters for frame 
(anholonomic) indices. 

We identify differential forms with covariant antisymmetric tensors. Given 
a pair of real covariant antisymmetric tensors P and Q of rank r we define 
their dot product a.s P ■ Q :— -^Pai a Qiii g"-^^^ . . .g°'''^'\ We also define 
\\Pr:=P-P. 

We define the action of the Hodge star on a rank r antisymmetric tensor R 
as {*R)a,.+i...a;i '■— (j'!)"^ ^"^ '"'^^ai...a3 where e is the totally antisymmetric 
quantity, £012 := +1- 

We use two-component complex-valued spinors (Weyl spinors) whose indices 
run through the values 1,2 or 1,2. Complex conjugation makes the undotted 
indices dotted and vice versa. 

We define the "metric spinor" 

ah ^db 



e., = 6,, = e- = 6-=(J qM (A.l) 



and choose Pauli matrices 



coab— ( 1 ) ~ ^'^ "^' (^-2) 
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fOl\ I fO-i\ 2 fl 0\ 3 ,. „, 

Here the first index enumerates rows and the second enumerates columns 



B. Nonlinear second order equations which reduce to pairs of linear 
first order equations 

Let n be an open subset of M" . We work with (infinitely) smooth vector func- 
tions il -> C™ writing these as columns of m complex scalars. In this appendix 
"vector" does not carry a differential geometric meaning because we are not in- 
terested in coordinate transformations. We use Cartesian coordinates x^, . . . , a:". 

Given a pair of vector functions u, w : J7 — > C™ we define their inner product 
in the standard Euclidean manner as (u, v) := J„ v*u dx^ . . . dx"' where the star * 
denotes Hermitian conjugation. This integral need not converge as we will be 
using it only for the purpose of defining the formal adjoint of a differential 
operator, see next paragraph. 

Let A± be a pair of formally self-adjoint (symmetric) first order linear partial 
differential operators (differential expressions) with smooth coefficients acting on 
smooth vector functions f2 -^ C™. We do not introduce any boundary conditions. 

Put 

L±(u) :=Re(uM±u). (B.l) 

It is easy to see that L±{u) is the Lagrangian density for the partial differential 
equation A±u = 0. Namely, if one writes down the action (variational functional) 
S±(u) := Jq L±{u) dx^ . . . dx" then the corresponding field equation (Euler- 
Lagrange equation) is A±u = 0. 

Let us now define a new Lagrangian density 

m .= /|(f^-^;\ (B.2) 

-L-i-(u) — L^yu) 

and corresponding action S{u) :— j' Liu) dx^ . . . dx". The field equation for the 
Lagrangian density (jB.2p is, of course, second order and nonlinear. 

Note that the notation in this appendix is self-contained and the Lagrangian 



densities (jB.ip . (jB.2p should not be confused with the Lagrangian densities (|4.5p . 
(|5.3p introduced in the main text (the latter have an extra subscript). 
The main result of this appendix is 

Lemma 1. Let u : f2 ^)- C" be a vector function satisfying the condition 

L+{u)^L^{u). (B.3) 

Then u is a solution of the field equation for the Lagrangian density L if and 
only if it is a solution of the equation A+m = or the equation A^u ^ 0. 

Proof. The explicit formula for the operator A± is 

A± = iB^dc, + '-{do^BD + C± (B.4) 
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where _B^ and C± are some smooth Hermitian m x m matrix functions and the 
index a runs through the values 1, . . . ,n. Substituting (IB.4|) into (IB.ip we get 

L±{u) = ^ [u*B'ldaU - {dcU*)Blu\ + u*C±u. (B.5) 

Now take an arbitrary smooth function /i : 12 — >■ M. Examination of formula 
(|B.5[) shows that 

L±{e'''u) = e2''L±(u). (B.6) 

We call the property (|B.6p scaling covariance. Scaling covariance is a remarkable 
feature of the Lagrangian density of a formally self-adjoint first order linear 
partial differential operator. 

Formulas (jB.2[) and (|B.6p imply that the Lagrangian density L also possesses 
the property of scalar covariance, i.e. L{e^u) = e'^^L{u) for any smooth h : fi ^ 
R. Thus, all three of our Lagrangian densities, L, L^ and L_, have this property. 

Observe now that if the vector function m is a solution of the field equation for 
some Lagrangian density C possessing the property of scaling covariance then 
£(u) — 0. Indeed, let us perform a scaling variation of our vector function 

u^ u + Su = u + hu = e''u + 0{h'^) (B.7) 

where /i : i7 — > R is an arbitrary "small" smooth function with compact support, 
h G C^(f2;R). Then = SjC{u) = 2 J hC{u) which holds for arbitrary h only 
iiC{u) =0. 

In the remainder of the proof the variation Su : f2 —> C™ of the vector function 
u : /2 — !> C™ is arbitrary and not necessarily of the scaling type (IB. 71) . The only 
assumption is that Su £ C§°{Q; C™). 

Suppose that u is a solution of the field equation for the Lagrangian density 
i-i-. [The case when u is a solution of the field equation for the Lagrangian 
density i_ is handled similarly] Then L+{u) = and, in view of formula ()B.3|) . 
L-{u) =/= 0. Varying u we get 



L^{u) — L^{u) J L+{u) — L^{u) 

6L+{u) = -5 I L+{u) 



S L{u) = ~S L+{u) . (B.. 



so 



We assumed that m is a solution of the field equation for the Lagrangian density 
L+ so SjL+{u) = and formula (jB.Sp implies that SjL{u) = 0. As the latter is 
true for an arbitrary variation of u this means that u is a solution of the field 
equation for the Lagrangian density L. 

Suppose that u is a solution of the field equation for the Lagrangian density 
L. Then L{u) = and formula (|B.2[) implies that either L^{u) = or L^{u) = 0; 
note that in view of ()B.3p we cannot have simultaneously i+(u) = and L-(u) = 
0. Assume for definiteness that L^{u) = 0. [The case when L-{u) = is handled 
similarly.] Varying u and repeating the argument from the previous paragraph 
we arrive at (jB.Sp . We assumed that w is a solution of the field equation for the 
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Lagrangian density L so SjL{u) — and formula (jB.Sp implies that SjL^{u) = 
0. As the latter is true for an arbitrary variation of u this means that u is a 
solution of the field equation for the Lagrangian density L^. D 

Remark 2. It may seem that the variational proof presented above is "insuffi- 
ciently rigorous". An alternative "completely rigorous" way of proving Lemma [1] 
is to write down the field equation for the Lagrangian density (jB.2p . (jB.sp ex- 
plicitly and analyze this second order nonlinear partial differential equation. The 
result, of course, remains the same, only the calculations become much longer. 

Remark 3. Examination of the proof of Lemma [1] shows that the fact that the 
differential operators A± are linear and first order is not important. What is 
important is that their Lagrangian densities possess the scaling covariance prop- 
erty (|B.6p . As the Lagrangian density (IB.2P possesses this property as well, our 
construction admits an obvious extension which gives a hierarchy of nonlinear 
partial differential equations which reduce to several separate equations. 

Example 1. Let us give an elementary example illustrating the use of Lemma [TJ 
Consider the pair of linear first order ordinary differential equations 

iu'±u^O (B.9) 

where u : M — )• C is a scalar function. Let us write down the corresponding La- 
grangian densities L±{u) = ^{uu' — uu') ± |up in accordance with formula (jB.ll) 

and form a new Lagrangian density —2L{u) = ( ""^7^" ) -|- \u\^ in accordance 

with formula (|B.2p . The latter gives the field equation (Euler-Lagrange equa- 
tion) 

' uu' — uu! \ [uu'Y — [uu'Y „ /„ ,„N 

- ' ' ^ ' ^ '—u + u = Q. (B.IO) 



2|u|2 J 4|u|4 

Lemma [1] tells us that a smooth nonvanishing function m is a solution of equa- 
tion (jB.lOp if and only if it is a solution of one of the two equations (|B.9p . Of 
course, this fact can be checked directly by switching to the polar representation 
u = re"*"^ where r : M ^ (0, +oo) and ^ : M ^ R. 
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